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Abstract 

In the ordered phase of 2D XY ferromagnet the dipole force in- 
duce a strong interaction between spin-waves in the long wavelength 
limit. This interaction leads to the transformation of the spin- wave 
excitation into a new soft-mode excitation in an intermediate range of 
wavelengths limited in magnitude and direction; and into an anoma- 
lous anisotropic diffusion mode excitation at long wavelengths. The 
dissipation of a spin-wave at short wavelengths is found to be highly 
anisotropic. 
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The conventional condensed matter theory deals with elementary exci- 
tations and their interactions. The excitations such as electrons, phonons, 
spin-waves etc. have a propagating, wave-like nature. The momentum p 
and the energy uo of a single excitation are related by the dispersion relation 
(spectrum): uo = e(p). A weak interaction changes slightly the spectrum and 
leads to a finite life-time of the excitations JT|. Effect of a strong interac- 
tion is not so universal. Migdal || has shown that a strong electron-phonon 
interaction renormalizes substantially electron velocities and, if it exceeds a 
critical value, destroys the electron quasi-excitations. Recently a strong in- 
teraction of electrons in 2 dimensions with gauge fields has been shown to 
result in non- Fermi-liquid behavior |J. 

In the long wavelength (hydrodynamic) limit quasi-excitations turn into 
classical modes, such as sound or spin-wave. Not only the propagating waves, 
but also particle, heat and spin diffusion can be considered as hydrodynamic 
modes. The interaction between these modes has been shown to be substan- 
tial in the critical region ||; in the dynamics of liquid crystals ||, f|; and 
in the dynamics of a Charge Density Wave interacting with impurities [^]. 
In all these systems the interaction leads to a drastic reconstruction of the 
dispersion relation. 

In this Letter we solve the experimentally feasible magnetic model, in 
which a strong interaction between spin-waves leads to the replacement of 
the propagating spin-wave by a diffusion mode and to appearance of a new 
soft-mode in a range of momentum. This model is the two-dimensional XY 
ferromagnet with the dipolar interactions between spins. 

A spin diffusion mode appears naturally in the paramagnetic phase and 
in the vicinity of the Curie point ||. We consider a low temperature ordered 
phase, where no diffusion is expected but rather a propagating and weakly 
dissipating spin-wave mode. In a 3D ferromagnet the exchange interaction 
between spin- waves vanishes in the long wavelength limit 0. The dipole 
force generates three spin-wave processes ^ [TD| and violates the total spin 



conservation law. This interaction is dominant in the dissipation of a spin- 
wave via decaying it into the two spin-waves or via merging it with other 
spin-wave. But dissipation is weak. 

In a 2D XY ferromagnet at low temperatures the dipolar interaction is 
relevant in the long wavelength limit, even despite of the low density of spin- 
waves. It was shown by one of the authors [|ll|] that dipolar force induces an 
anomalous anisotropic scaling of spin-spin correlations in the ordered phase. 
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In this article we find an analogous dynamical scaling. 

The dipolar force stabilizes the ferromagnetic long-range order [12], sup- 
pressing strong XY thermal fluctuation. Therefore, we represent the unit 
vector field of magnetization S by the two spin-wave fields - in-plane 0(x, t) 
and out-of-plane 7r(x, t) = S z : 

S = (— Vl — vr 2 sin0; y/l — ir 2 cos0; 7r) (1) 

where both n and <j) are small. The field ir(t, x) is canonically conjugated to 
the field <p(t, x). 

The Hamiltonian of a 2D XY ferromagnet contains three terms: the ex- 
change and the anisotropy energies: 

Hea = E( J P 2 |<M 2 + AK P | 2 )/2, 

LU p 

as well as the dipolar force term (see e.g. ]|iT| ): 



H = H E a + 9 ^ V-Aup + Vy (0 2 / 2 



w p 



up 
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/2|p|, (2) 



with the corresponding couplings being - the exchange constant J, the anisotropy 
A and the dipole constant g [13]; The Hamiltonian (Q) is written in terms 
of the Fourier transformed fields - <^ p and 7r wp , the abbreviation - (</> 2 /2) pw 
- denotes the Fourier transform of 2 (x, t)/2. The anisotropy is assumed 
to be weak: A C J. Therefore the spin S averaged over scales larger than 
J J/X turns into the plane. The bare spin-wave spectrum extracted from the 
quadratic part of the Hamiltonian (0) reads: 

e 2 (p)=c 2 (p 2 + PoP l/p), (3) 
where p = |p|, the characteristic wave- vector po = g/J and the spin- wave 



velocity c = v JX. In the region po ^ p ^ y X/ J spin-waves have the linear, 
acoustic-like spectrum. We call this spherical shell of momentum the ^4-shell. 

The non-linear part of the Hamiltonian (^]) Hint contains the three-leg 
and four-leg vertices. In particular the three-leg vertex is: 

3 

/(Pi, P2, P3) = Y.VixPiyl | Pi | (4) 
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Since the sum of momenta entering the vertex is zero it depends on two 
momenta. Further we denote it as /(pi,P2)- 

The dynamical equation associated with the Hamiltonian (□) reads: 



'J 2 e 2 (p) . uj \ , 8H int 

~i + Vvp- (5) 



X 




The thermal noise rj^p and the bare dissipation coefficient Tq 1 generate the 
stochastic dynamics in (|5]). The noise-noise corelation function is related to 
the bare dissipation coefficient: (|^ p | 2 ) = 2T/T . It vanishes at T = 0. We 
do not consider external sources of dissipation, spin-waves dissipate due to 
interaction generated by dipole force. Hence, we set T^ 1 = +0, and look for 
the generated dissipation: r _1 (c<;,p). 

The Green function G(u,p) is defined as the linear response of the mag- 
net to an external magnetic field with the same frequency and wave-vector. 
The l.-h.-side of Eq.[5] represents the inverse of the bare Green function 
A _1 Go 1 (a;, p). According to the Fluctuation-Dissipation Theorem, the imag- 
inary part of the Green function multiplied by the factor T/uj is the spin-spin 
correlation function D(u,p). For convenience we use the dissipation func- 
tion defined as: b(u,p) = X/T(u,p). The self-energy term E(a>,p) equals 
G , _1 (co',p) — G _1 (c<j,p) by definition. We notify the real and the imaginary 
part of the self-energy term as: £ = a 2 (u,p) — iu>b(uj,p). Thus, the Green 
function reads: 

G^^lj, p) = to 2 — e 2 (p) — a 2 (uj, p) + itob(uj, p), (6) 
whereas the spin-spin correlation function reads: 

D(u,p) = g (') 

[uj 2 — e 2 (p) — a 2 (u, p)] + u 2 b 2 (u, p) 

(we refer the factor T to vertices). We use the reduced temperature t = 
T/AttJ as a small and L = log(V JX/g) as a large parameters (the latter 
means large ^4.-shell). 

We employ the standard, so-called Janssen-De-Dominicis Functional method 



[|14|j to account for the non-linear terms in the stochastic Langevin equation 
(HD- This method generates a diagrammatic expansion in powers of the bare 
vertices, associated with H int . The main contribution to the self-energy is 
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given by the one- loop diagram shown in Fig. la, lb. Our theory is valid only 
if the temperature is small: 



t\og(VJX/g) = tL < 1. 
Under this condition the two-loop corrections are small and the diagram 



Fig. lb contributes to a negligible change of the spectrum (|3]) P3| . Such ne- 
glecting of the two-loop diagrams (vertex correction) was a major assumption 
in the so-called mode-coupling methods f| . Later we proove this assumption. 
Thus, the Dyson equation for our problem is as follows: 



E ( n>q ) = VAT/^/^/»(p,q) 



D(u,p)G(uj + n,p + q), (8) 

The functions b(u, p) and a(u, p) are even in both arguments 0. The imag- 
inary part of the self-energy is odd in frequency Q. Hence, the equation for 
the dissipation function reads: 

d 2 p r duj 



(27T) 

D(uj,p)D(uj + n,p + q). (9) 

The integrand in ([|) is positive. Thus, the main contribution to b(Q, q) 
comes from the region, where poles of the two D-functions coincide. The 
function D(u,p) has poles at: u ~ ±e(p), in the ^4-shell |TBJ. Following the 
terminology of the field theory we call the surface uj 2 = e 2 (p) the mass-shell. 
The dissipation in the ^4-shell is small and the D-function can be represented 
as a sum of 5-functions: 



± 



26 2 (p) 



where Au± — u ± e(p) measures the deviation from the mass-shell. After 
integrating u out from the Eq.|9| with the D-functions from ( p!Q|) we recover the 
Fermi Golden Rule for the probability of the spin-wave decay and scattering 
processes. 

Looking for the long wavelength quasi-excitations we need the self-energy 
at very small momenta q <^ p , which we denote as E . We expect quasi- 
excitations to be soft: Q <^ cq, and we restrict the wave-vector q of quasi- 
excitation to be directed along the magnetization: \q x \ <^ q. The essential 
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contribution to the integral in Eq|| comes from the internal momentum p 
being in the .4-shell and the internal frequency u = e(p). Taking the integral 
over u, with the D-function from Eq. (]TTJ|) , we find |T6[ : 

c 2 p\t f c A p 3 dp Q sin 2 (2ip) dip 



Aix J e 4 (p) Vt — cq cos ip + ib\ 

where b\ is the dissipation function b(u, p) of a spin-wave inside the ^4-shell 
and ip is the direction of the internal spin-wave: smip = p x /p. Note, that 
we added a f2-independent contribution from the diagram Fig.l.b to the self- 



energy [Ty. ReS vanishes in the static limit £7 = 0. 



If cq ^> b\ we make the integral over ip in Eq.|TT] to find: 

S (x) = c 2 p 2 tL cos 2 xexp(-2ix), (12) 
where x> defined by the equation cosx = £l/cq, measures the deviation 



from the mass shell. If q is so small that cq <^ b\ the Eq.|ll] implies the 
g-independent dissipation constant: 



c 2 p\tL J 



# sin 2 (2V>) 

4tt h(ip) ' 1 ' 



In this calculation we have used the fact that the dissipation of a spin-wave 
in the -4,-shell b\ depends only on the angle ip between the direction of mag- 
netization and the spin-wave wave-vector p, which we prove below. 

Anyway, we need bi(ip). An unusual feature of our theory is that the 
dissipation process in the ^4-shell is mediated by an off-mass-shell, virtual 
spin-wave. Indeed, the Eq.^| does not allow for decay or merging processes. 
Alternatively, the dissipation of a spin-wave in the v4-shell propagating in the 
direction specified with the angle ip (simp = q x /q) is mediated by an internal 
virtual spin- wave in (9), with a momentum p p q and a frequency 
uo < cp. The integration over uo with one of the D-functions in (|9|) taken in 
the form ([H]), leads to a following equation: 

= cW(0,q) ( rf2pD(e(p + q) _ e(q); p) (u) 
8rrq z J 

Since u = e(p + q) — e(q) we conclude that u = cp cos $, where $ = ip — ip is 
the angle between the vectors q and p. Invoking the definition of the angle 
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X for virtual spin-wave, we find that x — Now, we look onto Eq.|TJ] in 
more details: 



i. I i\ _ ^P& r f Sin2 1 COS t dpd(fi fltz\ 

b lW - ^ imj p2sin 2^ +po ^ 2 + SoW/c2 - W 



Note, that <p ~ p/po 1 hi (P^D- Thus, x = V*- m other words, the 
dissipation of a short wavelength spin-wave propagating in the direction tp is 
determined by the scattering on the long wavelength virtual spin-wave, which 
lies on the specific distance off the mass-shell: uj/cp = cosip. The integration 
over p in fll5|) is confined towards the crossover region: p ~ p c = p^yftL. 

Substituting £o(V0 from Eq.|TJ into Eq.[L5| we find the anisotropic dissi- 
pation of a spin-wave mode in the ^4-shell: 

, / ,x n im sin 3/2 (2^) sin(W2) , , 

L 1 / 4 cos?/; 

where the direction of the spin- wave is limited to the fundamental quadrant: 
< ijj < tt/2. We found: f3 l = T 2 (1/A)/AV2k » 1.31. 



Let us return to very low momenta p <C b\/c. Plugging Eq.|l6| into EqJJJ 
one finds: &o — f3ocpot 1 ^L 5 ^ 4 , where (3q rs 1.24 was found numerically. The 
condition cpou ~ b\ defines the crossover wave-vector: p DM ~ p t 3 ^ / L 1 ^ 4 , 
between the self-energies Eq.|12] and Eq.|l^. 

The dissipation functions (|l^Jl~3| and [16]) conclude the self-consistent so- 
lution of the Dyson equation (||, §)• 

We verify that the two- loop correction (see Fig.lc) is negligible. Note, 
that the main contribution to the diagram Fig.lc comes, if the two internal 
momenta pi and p 2 are restricted to the ^4-shell. Inside the ^4-shell the Green 
and the D-functions live on the mass-shell. But, the prohibition of the three 
spin- wave processes leaves the Green function G{lo\ — u 2 , Pi — P2) to be off 
the mass-shell. The Green function off the mass-shell is small in temperature. 
A simple counting shows, then, that the two-loop dissipation function is 
6' = 6 (l+t 1/4 O(t)). The funct ion b'<y(ip,q) — bi(ip), which represents the 
two-loop corrections for bi, is small in t 1//4 , and is also small in the ratio 
Po/q- 

Having the explicit expression for the self-energy we can analyse the dis- 
persion relation - u 2 = e 2 (p) + £(u;, p) - in the range of small u and p. New 



results are expected for the region p < p c = p^ytL in which S becomes 
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comparable with e 2 (p). In a range of momentum Pdm <C p <C p c and angles 



ip <C ypotL/p we find a new propagating soft mode with the dispersion: 

u = cp{p 2 + p pip 2 ) 1/2 /p VtL (17) 
The dissipation of the soft mode grows to the boundary of the region and 



becomes of the order of its energy at ip ~ ^potL/p or p ~ Pdm- There is 
no soft mode beyond the indicated range. The spin-wave mode persists at 
p > poVtL. In a range p <C Pdm and small angles a new diffusion mode 
occurs with the dispersion: 



UJ 



-ie^r^L^/Pocpa (18) 



The angular range of the diffusion mode increases with decreasing p and 
captures the entire circle at p < PotL. 

At even smaller wavelengths p < Pa ^ Pdm the interaction between dif- 
fusion modes should be taken into account ( pa is the anomalous diffusion 
set up wave- vector [H]). The problem can be solved by the renormalization 



group method [pT |. The growing interaction, although leaving invariant the 
diffusive nature of the spin propagation, changes the dispersion. For the prop- 
agation along the spontaneous magnetization it is: uj oc —ip A7 ^ 27 ) whereas for 
the propagation in the perpendicular direction it is: to oc —ip 47 ^ 36 . This is a 



dynamic analog of the non-Gaussian fixed point found in |T . 

In conclusion, we discuss the experimental feasibility of the above consid- 
ered effects in the epitaxial magnetic films. The main difficulty is that all of 
them are confined to rather long wavelengths. Therefore even weak in-plane 
anisotropy can suppress them. One need to use substrates with the six-fold 
symmetry axis. The six-fold anisotropy is much weaker than the tetragonal 



one. Moreover, it was predicted in 18 that there exists a temperature inter- 



val in which the hexagonal anisotropy vanishes at large distances. A proper 
substrate is, for example, a (111) face of Cu, Au etc. Unfortunately the 
epitaxial growth of ferromagnetic films on this faces has not been realized so 
far. In a recent work ]TJ| a growth of a ultrathin Ru/C(1000) film has been 
reported. It is the in-plane ferromagnet. We hope that further sophistication 
of the experimental technics will allow to observe the predicted dynamical 
behavior. 

We are indebted to M.V. Feigelman, E.I. Kats and V.V. Lebedev for 
useful discussions and indicating some references. 
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Figure Caption 

• FIG.l. a) Main one-loop self-energy diagram; b) One-loop four-leg 
vertex diagram; c) Two-loop self-energy diagram. Momenta of internal 
lines are indicated. 
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